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We show that electrically and magnetically frustrated Josephson junction arrays (JJAs) realize 
topological order with a non-trivial ground state degeneracy on manifolds with non-trivial topology. 
The low-energy theory has the same gauge dynamics of the unfrustrated JJAs but for different, 
"fractional" degrees of freedom, a principle reminescent of Jain's composite electrons in the fractional 
quantum Hall effect. 
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The idea of quantum order has its origin in the discov- 
ery of the fractional quantum Hall liquids [l|, that are 
incompressible and exist only at some "magical" filling 
fractions. The internal order characterizing these states 
is a new type of order, called topological order, different 
from any other previously known type of order and it is 
not based on spontaneous symmetry breaking. Quantum 
order describes in general the zero-temperature proper- 
ties of a quantum ground state and characterizes uni- 
versality classes of quantum states, described by com- 
plex ground state wave-functions [2j. Topological order 
is a special type of quantum order whose hallmarks are 
the presence of a gap for all excitations (incompressibil- 
ity) and the degeneracy of the ground state on manifolds 
with non-trivial topology Q • In the case of the fractional 
quantum Hall effect, which is a (2+l)-dimensional sys- 
tem, another hallmark is the presence of excitations with 
fractional spin and statistics, called anyons @. The long- 
distance properties of these topological fluids can be ex- 
plained by an infinite-dimensional W^i+oo dynamical sym- 
metry 0], and are described by effective Chern-Simons 
field theories [j| with compact gauge group, which break 
P-and T-invariance. In fact, the simplest example of 
a topological fluid Q is a ground-state described by a 
low-energy effective action given solely by the topological 
Chern-Simons term S = k/4n f d 3 x A^ ua d u A a for 
a compact U{\) gauge field whose dual field strength 
— e^ va d v A a is interpreted as the conserved mat- 
ter current. In this case the degeneracy of the ground 
state on a manifold of genus g will be k 9 (or (£;ifc2) 9 
if k = k\/k2 is a rational number). Topological order 
in conventional superconductors has been investigated in 

In Q we proposed a new superconductivity mechanism 
which is based on a topologically ordered ground state 
rather than on the usual Landau mechanism of spon- 
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tancous symmetry breaking. Contrary to anyon super- 
conductivity it works in any dimension and it preserves 
P-and T-invariance. The topologically ordered supercon- 
ductors we proposed have a long-distance hydrodynamic 
action which can be entirely formulated in terms of gen- 
eralized compact gauge fields, the dominant term being 
the topological BF action [10|. In (2+1) dimensions, it 
reduces to a mixed Chern-Simons model [Til ] . 

A first concrete example of topological superconduc- 
tors is provided by JJAs. In [l2j , we have proven that pla- 
nar JJAs can be exactly mapped onto an Abelian gauge 
theory with a mixed Chern-Simons term (BF-model) 
which complete determines the low energy dynamics. 
The Abelian gauge theory exactly reproduces the phase 
diagram of JJAs and the insulator/superconductor quan- 
tum phase transition at T = [13[. JJAs, and fully 
frustrated J J A ladders |l4j . have also been recently 
considered by several other authors [HI, [lj| as devices ex- 
hibiting topological order. For planar unfrustrated JJAs, 
however, the topological fluid is described by a k = 1 
mixed Chern-Simons model and, thus, there is no degen- 
eracy of the ground state [1, [12] • 

In [13] we argued that frustrated JJAs may support a 
topologically ordered ground state with non-trivial de- 
generacy on higher genus surfaces. In this paper we 
give an exact derivation of the gauge theory that de- 
scribes Josephson junction arrays in presence of electric 
and magnetic frustration. We will then show, for ratio- 
nal frustration, how this can be rewritten as the same 
dynamics of the unfrustrated model but for different de- 
grees of freedom, a principle analogous to Jain's compos- 
ite electrons familiar from the fractional quantum Hall 
effect. The only difference is the Chern-Simons coeffi- 
cient, determined by the denominator of the frustration: 
as we anticipated this leads to a ground state degeneracy 
for arrays with non-trivial topology. 

We shall consider JJAs fabricated on a square planar 
lattice of spacing I = 1 made of superconducting islands 
with nearest neighbours Josephson couplings of strength 
Ej [TU ■ Each island has a capacitance Co to the ground; 
moreover, there are nearest neighbours capacitances C . 
To implement a torus topology we impose doubly pe- 
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riodic conditions at the boundary of the square lattice. 
The Hamiltonian describcing the system is, in the limit 
C>C : 



comes: 



H 



J2 N 2 E c pa^Po 



Ej (1 - cos N (Ai$)) 



where Eq = q^/2C. The phases $ are quantum- 
mechanically conjugated to the charges Q on the islands: 
these are quantized in integer multiples of N (Cooper 
pairs for N = 2): Q = q e Np , po 6 Z where q e is the 
electron charge. In [l2j we have shown that the zero- 
temperature partition function of JJAs may be written 
in terms of two effective gauge fields (vector) and 
(pseudovector). In the low energy limit the partition 
function is 



= J2 ( SM o) J VA ^J VB v exp(-S) , 



{Qo> 
{Mo) 



s 



-- J dtj^-i^ A 1l K ia) B v + 
-iAoQo + iB Q M . 



(2) 



This form of the partition function holds true also with 
toroidal boundary conditions. The gauge fields, with 
compactification radius 2-7r, embody the original degrees 
of freedom through their dual field strengths p^ oc K^ U B U 
and g M oc K^ V A V representing, respectively, the con- 
served charge (vector) current and the conserved vortex 
(axial) current. is the lattice Chern-Simons term 

[ll], defined by K 00 = 0, K 0l = — eydj, K i0 = S l e l jd j 
and Kij = — S^e^do, in terms of forward (backward) shift 
and difference operators Si (Si) and di (di). Its conjugate 
K^ v is defined by K Q0 = 0, K 0l = S,t ,,d,. K M = e lj d J 
and = —SjCijdo. The two Chern-Simons kernels 
and K^ lL , are interchanged upon integration (summation) 
by parts on the lattice. 

The topological excitations are described by the 
integer-valued fields Qo and Mq and represent unit 
charges and vortices rendering the gauge field compo- 
nents Aq and Bq integers via the Poisson summation 
formula; their fluctuations determine the phase diagram 
[l2| . In the classical limit the magnetic excitations are di- 
lute and the charge excitations condense, rendering the 
system a superconductor: vortex confinement amounts 
here to the Meissner effect. In the quantum limit, the 
magnetic excitations condense while the charged ones 
become dilute: the system exhibits insulating behavior 
due to vortex superconductivity accompanied by a charge 
Meissner effect. 

In presence of an homogeneous charge frustration, rep- 
resented by external charges pf on the islands fl} be- 



H= N2E c (Po ~ Pf) (PO - Pf) 

X 

^2 Ej (l -cos N(Ai$)) , 



(3) 



Following the same steps as in [12j we can write the par- 
tition function for the JJAs as: 



Z=Yj f Va nf V K exp(-S) , 
{Qo} J 

S = [ dt V -i2n a„K„ u b v + + 
+N 2 l a E c (po-Pf) — ^ (po-p f ) +i2Tra Q Q , (4) 

with p^ ee K^by , & M 6 R , o M = K^a^ , a M e Z 

In this representation K^b v represents the conserved 
three-current of charges, while K^a^ represents the con- 
served three-current of vortices. a M and 6 M are compact 
gauge fields with period 1. The third term in the action 
([JJ contains two parts: the longitudinal part (pf) de- 
scribes the Josephson currents and represents a kinetic 
term for the charges; the transverse part (pf) can be 
rewritten as a Coulomb interaction term for the vortex 
density qo by solving the Gauss law enforced by the La- 
grange multiplier bo. 

The partition function ^ displays a high degree of 
symmetry between the charge and the vortex degrees of 
freedom. The only term which breaks this symmetry 
(apart from the integers Qo) is encoded in the kinetic 
term for the charges (Josephson currents). Following 
[l2| we introduce the self-dual approximation of JJAs 
by adding to the action in ([4]) a bare kinetic term for the 
vortices Yl x N ^ E qf The coefficient is chosen so that 
the transverse part of this term reproduces exactly the 
Coulomb term for the charges upon solving the Gauss 
law enforced by the Lagrange multiplier ao- 

In order to exactly reproduce the Coulomb term, in- 
cluding the frustration, we need to add an extra ilixaoPf ■ 
The action then becomes: 



Ssd = J dt ^2 -* 27r a^K^bv 



P, 



-\ 2 



2l Ej m E a 
i2iraoPf + i2iraoQo + i2nboMo , (5) 

where we are also forced to introduce new integers Mq via 
the Poisson formula to guarantee that the charge current 
Kfjtvbv remains an integer (note that, without the kinetic 
term for the vortices, both conserved currents are inte- 
gers: indeed, the summation over {Qo} makes ao , and 
then the summation over {ao} makes K^ v b v an integer). 
Since the properties of the system are periodic with re- 
spect to the charges, we can restrict —1/2 < pf < 1/2 
|19j |: henceforth we have Pf = p/q. Since we are inter- 
sted only on the ground state, we will consider only to 
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the low energy limit. In this limit the partition function 
becomes: 



Zle 
Sle 



= ^2 S {Mo) I Pa„ j X>6 M exp {-Sle) , 

(Co) 

— / dty^ — i2tt a^Kpybv 

J X 

P 

+i2naoQo + tlith^M^ + z27rao — . (6) 



The effect of charge frustration, the last term in |6]), 
can be reabsorbed as follows. We first change the com- 
pactification radius of the gauge field ao from the stan- 
dard 1 to q. As usual, this corresponds to admitting 
fractional charges N/q in the model (in order to keep 
with standard notation we rescale a M — * qa^ , where a M 
has again period 1). Correspondingly we have to rescale 
the magnetic topological excitations Mq — * qMo in or- 
der to be able to absorb integer shifts in the rescaled 
a^. This corresponds to having a condensate of fraction- 



ally charged particles N/ q while the magnetic excitations 
have fractional (in units of the magnetic flux quantum 
2nq/N) vortices 2ir/N. Since the magnetic topological 
excitations qMo represents tunneling between vortex sec- 
tors differing by q units 2?t/N we have a Z q theory in 
which fractional vortices and charges have mutual frac- 
tional statistics 2Tr/q [2(|. For a = 2 such mutual anyons 
where considered by Kitaev [21| in his model of topolog- 
ical quantum computation. At this point the remaning 
frustration effect ca be easily accounted by a redefinition 



6 = b ;h = h + Af* 



where we have written, without loss of generality, the 
frustration term pf as pf — K^Af^ . The action ([6]) 
becomes 

Sle = / dty^— i2nq a^K^yby 

+i2TTqa Q + i2nqb M + i2Trqa K 0i A^ t . (7) 



obtain the low energy partition function: 
Zle = <H M °) JvA^JvB^ exp (-Sle) 



{Qo> 
{Mo} 



Sle = -i J dt^-^-A^KnyBu 
iqB M + iqAoQ . 



(9) 



This is the same action as for the unfrustrated case (same 
periodicity of the gauge fields by construction) with the 
only difference of a Chern-Simons coefficient q deter- 
mined by the denominator of the frustration. In analogy 
to the Jain construction familiar from the QHE we have 
the same dynamics but for different degrees of freedom. 
This, however, affects the ground state degeneracy that 
is given [13] by (k) 9 = q on the torus. 

The magnetic frustration can be treated in a similar 
way. In presence of a uniform magnetic frustration the 
action Q] becomes: 



H 



with AS 



J2 N 2 E C po-^Po + 

X 

]T £./(!- cos (NA^-A^)) , (10) 



time independent and such that 



Epiaquette ^f* = 27T f where 2n f is the flux piercing the 
elementary plaquette. Due to the periodicity properties 
of the system [19|, we can restrict 0<p/<l/2 and 
rewrite: / = p/q. Following again the same steps as in 
1121. we arrive at 



Z 



s 



= J dt^-iNp^ (a^ + ^v, 



i Pt A° 



9 1 

+iN 2 E c l Po—rPo + 



Pi 



-A" u 2l Ej 



(11) 



Since Af is time- independent we have ciiKijbj 
ciiKijbj and we get 



ale = 8 (Mo) / Va^ / Vb^ exp {-Sle) , 

Wo) 

{M } 

Sle = -i I dt*^22Trqa K 0l b i + 2irqa l K lQ b + 



2irqa,iKijbj + i2irqboMo + i2irqcioQo . (8) 



Finally, by rescaling A^ = 27ra M , = 2-Kb^ in ©, we 



Introducing the compact gauge fields and 6 M , both 
with period 1, such that K^ v b u represents the conserved 
three-current of charges, and K^a^ represents the con- 
served three-current of vortices, we can rewrite the action 
in (jnH as 

+ib K oi Af ct + i2Tra Q + i2TTb M , (12) 

where we have introduced the two integer fields Mq and 
Qo as in ©. 
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The low energy partition function is: 



J LE 



s 



LE 



= HM ) fva^ jvb^ exp {-S LE ) , 

{Qo> 

= / dt ) -i2-K a^K^vby 

J X 

+i27raoQo + i27r&o-A^o + *27r6 — ■ (13) 



Following the same steps as for charge frustration we 
reabsorb the effect of magnetic frustration, the last term 
in (|13p. as follows. We first change the compactifica- 
tion radius of the gauge field bo from the standard 1 to q; 
this corresponds to admitting fractional fluxes 2ir/Nq.Wc 
then rescale b^ — > qb^ , where b^ has again period 1. 
Correspondingly, we have to rescale the charge topolog- 
ical excitations Qo — * qQo hi order to be able to absorb 
integer shifts in the rescaled 6 M . In this case we have frac- 
tional magnetic fluxes 2ir/Nq while the charge excitations 
have fractional (in units of the electric charge quantum 
Nq) charge N. Since the charge topological excitations 
qQo represents tunneling between charge sectors differ- 
ing by q units N we have a Z q theory, in agreement with 
the result found in 23], in which fractional vortices and 
charges have mutual fractional statistics 27r jq. We now 
change variables: 

a = a ; a i = a,; + A^ xt , 

with A| xt time-independent and rewrite the action in JT3]) 
as: 

Sle = dt^^ -iZnq a^K^K + i2qTra Q 

J X 

+i2qirb M . (14) 



By rescaling = 2-kcl^, B p 
the low energy partition function: 



2-irb^ in (1141) we obtain 



Zle = <K M °) [ VA v t VB » ex P (~ S le) , 

{Qoi 

Sle = -ijdt^2 t^A^K^B v + iqB M + iqA Q$) 



Again this is the same action as for the unfrustrated case 
(same periodicity of the gauge fields by construction). 
The only change is again the Chern-Simons coefficient q 
determined by the denominator of the frustration, giving 
a ground state degeneracy (k) 9 = q on the torus [221 ]. 

Planar frustrated Josephson junction arrays 
with toroidal boundary conditions represent the 
simplest example of a superconductor that exibits 
topological order with non-trivial ground state 
degeneracies. Furthermore they are accessible ex- 
perimentally [25] and thus they may be used to 
test the topological superconductivity mechanism 
proposed in [171 ]. 

The result we have obtained is very similar to the one 
obtained in [24| for the the frustrated Heisenberg model. 
Here the frustration is obtained by introducing next to 
nearest-neighbors interactions. For the chiral spin phase, 
the low energy limit is described by a Chern-Simons the- 
ory with coefficient k = q if the effective flux threading 
the plaquette is <j> = 2tt^. This is true also for (j> = it, 
because due to the frustration there is an effective flux 
tt/2 on triangular plaquettes. This hints to a possible 
equivalence between geometric frustration and frustra- 
tion induced by an external field 
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